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We have studied a one dimensional Hubbard superlattice with different Coulomb correlations at 
alternating sites for a half-filled band. Mean field calculations based on the Hartree-Fock approxi- 
mation together with a real space renormalization group technique were used to study the ground 
state of the system. The phase diagrams obtained in these approaches agree with each other from 
the weak to the intermediate coupling regime. The mean field results show very quick convergence 
with system size. The renormalization group results indicate a spatial modulation of local moments 
that was identified in some previous work. Also we have studied the magnetoconductance of such 
superlattices which reveals several interesting points. 
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I. INTRODUCTION 

The study of electronic correlation remains at the focus 
of recent theoretical interest due to the development of 
different novel materials of which the metallic multilayers 
[l[ form a very important component. The oscillation of 
exchange coupling between magnetic layers Q and the 
appearance of giant magnetoresistance [3| are among the 
exciting features of the magnetic multilayers, e.g. the 
layered Fe/Cr structures. An intense theoretical attempt 
has been made to understand the magnetic behavior of 
such systems [1, d, H, 0, Hi • In these pioneering works in 
this field a simple generalization of the one dimensional 
Hubbard model was proposed to investigate the role of 
electronic correlation in one dimensional superlattices. 
This model consists of a periodic arrangement of Nu sites 
in which the on-site Coulomb correlation {U) is repulsive 
followed by A'o sites with no on-site interaction [U = 
0). It was found that such a model gives rise to some 
interesting features, in sharp contrast with the magnetic 
behavior observed in an otherwise homogeneous system 
(conventional one band Hubbard model) [J]- 

However, most of the previous works in this model em- 
ployed exact diagonalization of finite systems typically 
having 8 — 24 sites [3, Q. Some of these works reveal 
preferential distribution of local moments on sublattices 
and observe the suppression of spin density wave (SDW) 
order 3. Some others deal with problems like metal- 
insulator transition (MIT) |5| or formation of charge den- 
sity wave in such systems. Density matrix renormal- 
ization group (DMRG) calculations are also performed 
to overcome limitations on system size imposed by the 
technique of exact diagonalization; and, these consider 
system sizes ^ 48 — 150 Q- Very recently a generaliza- 
tion of the said model has been considered where, instead 
of two types of sites with on-site correlation parameters 
U > Q and /7 = respectively, one takes into account 



two different values of [/(> 0) at adjacent sites 0; how- 
ever, a detailed study of this generalized model is yet to 
be worked out. In view of the wider applicability of this 
model to diverse experimental systems, we consider here 
a preliminary study of this alternating Hubbard model 
(AHM) in one dimension. Also we aim at observing the 
behavior of the same when the system size is reasonably 
large. 

In the present study, we investigate the ground state 
properties of the AHM in one dimension and for a half- 
filled band by using a Hartree-Fock approximation (HFA) 
together with a real space renormalization group (RG) 
calculation. These two techniques, complemented by one 
another, was found to be very successful in studying sim- 
ilar cases in recent past [10]. Apart from constructing 
the ground state phase diagram we also investigate the 
magnetoconductance of a finite chain within the HFA 
scenario. In sec. II we introduce the model and give the 
HFA calculations. Sec. HI describes the results obtained 
in HFA. Sec. IV contains some details of the RG scheme 
while sec.V shows the RG results together with a com- 
parison between the same and that obtained in the HFA. 
In sec. VI we present the results on magnetoconductance 
of the model and sec. VII summarizes the present work. 



II. THE MODEL AND THE HARTREE-FOCK 
CALCULATIONS 

Our model is defined on a one dimensional Hubbard 
chain of N (even integer) sites consisting of two sublat- 
tices. The model Hamiltonian is, 

i,(T 

+UA^ni^'jni^l + UB^ni^^ni^l (1) 
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The two sublattices constructed out of odd and even 
numbered sites are labeled by A and B respectively. 
cj^(c,,(j) is the creation(annihilation) operator for an 

electron with spin a at the i-th site. Ui^a = c\^^Ci^a, 
and rii = Ui^cr is the number operator at the i-th 
site, t is the hopping integral between nearest neighbor 
sites. Ua, Ub are the on-site Coulomb repulsion energies 
on the sites corresponding to two sublattices A and B 
respectively. 

We decouple the Hamiltonian within the HFA, which 
is expected to work at least in the weak coupling regime. 
We define two parameters, number of electron TV, and 
the magnetization Mj at the i-th site, where. 



(2) 



These lead to a decoupled Hamiltonian: 

i,(T 
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--AUA + UB){Nf-Mf) 



Mi)ni,t)] 



(3) 



Now the Hamiltonian can be divided into two parts for 
two types of spins, i.e. H = H-^ + H^. In an unrestricted 

HFA one diagonalizes and in a self-consistent 
manner to obtain the single particle energy levels. The 
ground state can be constructed by filling up the energy 
levels from both the up and the down bands upto the 
Fermi level. 

One can define the spin and the charge density order 
parameters, c and s respectively, by 



-y 



-n, 



We consider a half-filled chain with periodic boundary 
condition. One can easily check by using an unrestricted 
HFA calculation that all sites corresponding to a given 
sublattice become equivalent under a periodic boundary 
condition. This leads to a simplification of the formulas 
for c and ,s given above. We, therefore, use these simpli- 
fied forms of the charge density order parameter c and the 
spin density order parameter or antiferromagnetic order 
parameter s as given by. 



s = ^(nsj -ns.i -»^A,T (4) 

It is to be noted here that |c| = 1 for a perfect "chess- 
board" -type CDW and \s\ = 1 for a perfect Neel-type 
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FIG. 1: Plot of the order parameters c and s as functions of 
Us/t for (a) Ua = 0.4, (b) Ua = 1.2, and (c) Ua = 2.5 (scale 
of energy is chosen by setting t = 1.0) for A'' = 100. The 
dotted line corresponds to the charge order parameter c while 
the solid line shows the spin order parameter s as calculated 
from HFA. 



antiferromagnetic SDW state. We now investigate the 
dependence of these order parameters on the values of 
Ua and Ub- 



III. THE RESULTS OF HARTREE-FOCK 
CALCULATIONS 

In fig.l we plot the order parameters c and s as func- 
tions of UB/t for different fixed values of UA/t (for a 
chain having sites N = 100). Keeping UA/t at a fixed 
value (> 0) and varying Us/t from to higher values, 
we find that initially the system is charge ordered. The 
electrons tend to localize at the sites with lower Coulomb 
repulsion energies, keeping the other sites vacant. As a 
result a charge density wave (CDW) is formed. In this 
regime, the charge density order parameter c assumes a 
high value while the spin density order parameter s is 
zero. As we gradually increase Us/t keeping UA/t at 
the same fixed value, we see that there is a gradual fall 
in the value of c. At a certain value of f/s/i there oc- 
curs a sharp rise in s which now takes over the value of 
c. For small values of UA/t{« 1) the transition occurs 
at the point of homogeneity, i.e. at Ub = Ua- In this 
case, however, a further increase in UB/t suppresses the 
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UbA 

FIG. 2: Plot of the energy gap /S.hfa, as found from HFA, as 
a function of Us/t for A'' = 300. We have used the dotted Une 
for Ua = 0.4, the sohd Une for Ua = 1.2, and the dashed Une 
for Ua ~ 2.5 (scale of energy is chosen by setting t = 1.0). 
Note that the zeroes in the gap match with the transition 
points identiUed from fig. 1. 

spin order again, and the CDW sets in. Therefore, the 
SDW is found to form only at a singular point Ua ~ Ub 
which is in agreement with the known result for the "ho- 
mogeneous" limit [llj. The situation becomes different 
for larger values of UA/t. For an intermediate value of 
UA/t{^ 1), we find two transitions: one from a CDW to 
an SDW and then from the SDW to a CDW again; this 
can be identified by two successive cross-overs in the c- 
and the s-curves (fig. lb). It is interesting to note that the 
charge order vanishes only at the point of homogeneity. 
For large UA/t{>> 1) there appears only one transition 
from a CDW state to an SDW state at a specific value 
of UB/t{< UA/t). Here also, the charge order param- 
eter vanishes at Ua = C^b, and then rises slowly with 
UB/t{> UA/t). However, in this region, the spin or- 
der parameter always dominates over the charge order 
parameter (fig.lc). The phase transitions occurring at 
points of cross-overs of c and s can further be explored 
by studying the gap in the spectrum at the Fermi level. 
The energy gap Ah fa at the Fermi level of a system 
containing n electrons can be estimated from 

Ah FA = En+l — En 

where, E„ is the ground state energy of (3) for a system 
of n particles. In fig. 2 we plot the energy gap Ah fa 
as a function of Un/t for different fixed values of UA/t 
at half-filling. Dips in the curves of Ahfa match with 
the corresponding values of Ub /t that were identified as 
points of phase transitions in fig.l. It is to be noted 
here that for large values of UA/t the energy gap sharply 
increases after the CDW/SDW transition has occurred 
(fig. 2). In this regime the slow increase of the charge 
order parameter does not modify this behavior. 

Identification of the points of phase transitions by the 




D 
W 




2. 4. 6. 8. 10. 

UaA 



FIG. 3: The HFA phase diagram of the AHM for different 
system sizes (the dashed line for N — 100, and the solid 
line for N = 200) in the UA/t - Ub /t plane. The plot is 
symmetrical about the "line of homogeneity" Ua = Ub which 
is shown by the dotted line. 



methods mentioned above enables us to draw the phase 
diagram of the model (1) on the UA/t—UB/t plane within 
the HFA. 

In fig. 3 we plot the phase diagrams for different N val- 
ues to note the quick convergence of the present mean 
field results with system size. For small values of Ua and 
Ub the antiferromagnetic phase actually occurs along 
the line of "homogeneity" {Ua = Ub). For higher val- 
ues of Ua and Ub we obtain a centrally located broad 
SDW region together with two charge ordered phases lo- 
cated near the axes. The phase diagram turns out to be 
perfectly symmetric about the line of homogeneity along 
which the system is already known to be antiferromag- 
netic from exact calculations [ll[. In the limit Ub >> Ua 
{Ua » Ub) the CDW/SDW transition line above (be- 
low) the hne of homogeneity bends towards the UB/t 
{UA/t) axis. 



IV. THE RENORMALIZATION GROUP 
CALCULATIONS 

Next we appl y a real space Renormalization Group 
(RG) technique [ij to the same model keeping in view 
the success of this technique to similar 1-d systems of 
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correlated electrons [l^iliS, [3 recent past. It is to be 
noted here that in order to achieve a closed parameter 
space under the present scheme of RG iteration, we must 
generalize the model (1) as follows: 



(5) 
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where, ca (es) refers to the site energy of a site belonging 
to the A (B) sublattice. We start with = es = which 
makes (5) equivalent to (1). However, in course of the RG 
iteration non-zero values of eA and es may appear. 

For implementing the RG transformation the whole 
chain is now divided into cells containing three sites each. 
Since our system is a bipartite lattice with two types of 
sites A and B, there will appear two types of cells ABA 
and BAB. At the renormalized length scale, we identify 
the ABA cells as the new (renormalized) A-type sites and 
the BAB cells as the new i?-type sites. There are four 
different 'on-site' states for each site |0), |-|-), |— ), | H — ). 
We diagonalize the cell-Hamiltonian (for both types of 
cells) , and among the eigenstates of the cell-Hamiltonian 
we retain only four low-lying states at each iteration, for 
construction of the RG recursion relations. We are inter- 
ested in the half-filled ground state; so we retain the low- 
est energy states in the subspaces {n = 2, S = Sz = 0}, 



{n — 3,S 
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±i} and {n = 4,5 



Sz 



0} 



of each type of cells. Here n, S, and Sz denote the total 
number of electrons, the total spin, and the z-component 
of the total spin respectively. These four states in a cell 
can now be identified as the 'renormalized' on-site states 
|0')> !+')> |-') and I + -') respectively. 

To find the renormalized hopping matrix element, the 
matrix elements of c^(A) and c^(i3) between 'renormal- 
ized on-site states' are calculated, where is the annihi- 
lation operator of the electron with spin a at the bound- 
ary site of the cell and A or _B in the parenthesis denotes 
the type of the cell. Let for an ABA type cell 

(0'|c?(A)|+') = Ai(A) , 

(-'|c?(A)| + -') = A2(A) , 

and for a BAB type cell 

(0'\c\{B)\+'} = X,iB) , 

(-'14(5)1 + -') = A2(i?) . 



(6) 



(7) 



Our system possesses spin reversal symmetry, so the ma- 
trix elements for cj 's will be as same as that of C| 's (ex- 
cept for a fermionic sign change in the value of A2). But 
due to the lack of particle-hole symmetry, Xi{A) ^ X2{A), 
and Xi{B) 7^ A^fS). At this stage we introduce an ap- 
proximation [13, [13, [iH by defining. 

A(^) = VAi(A)A2(A) 



which leads to 



X{B) = ^Xi{B)X2{B) 



c^(r) = A(r)c;(r) 



(8) 



(9) 



where T ~ A,B and a —'\, |. So the effective renormal- 
ized hopping becomes 



t' = X{A)X{B)t 



(10) 



The intra-cell Hamiltonian, restricted to the subspace 
of the four states |0'), |-|-'), |— ') and | H — '), can now be 
written in terms of the new cell-fermion operators [l^ as 

H' = Eo, + {E+, - E„,)in\+n[) 

+ {E+^. + Ea, - 2E+,)n\n\ (11) 

where E^ /, E+i and Eq' are the lowest energies of the 

subspaces corresponding to 4, 3 and 2 particles respec- 
tively. From this we can easily identify the renormalized 
on-site quantities 



and 



Ui. = E+^,{T)+Eo'{r)-2E+,{r) (12) 



e'r^E+,{r)-Eo'{T) (13) 



where T = A,B. 

The ground state energy per site is computed from the 
sum 



En = 



^ 00 



l4"'('l) + 4"'(B)l 



n=l 



3" 



(14) 



where n denotes n th stage of iteration. We also calculate 
the local moment Lq defined by 



{n 



In absence of particle-hole symmetry this leads to a re- 
cursion relation of the form flO'l 



Ln 



a + bL'o + cP' 



(15) 



where, P = [ni + n^) (n| + — 1) obeys a similar re- 
cursion relation given by 



P = d + eL'„ + fP' 



(16) 



with a,b,c,d,e, and / are obtained from the matrix el- 
ements of Lq and P between the truncated basis of the 
cell Hamiltonian. The operators Lq and P are considered 
for the central site of the cell to minimize the boundary 
effects [lol [T2j. However, in implementing the recursion 
relation it is to be noted that a recursion for Lq at a i? 
type site (being at the middle of an ABA cell) will in- 
volve Lq and P' pertaining to a renormalized A type site, 
because under the RG transformation the ABA cell a 
'renormalized' A site. Similar consideration arises for P 
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FIG. 4: RG flow diagram in an effective parameter space 
Us/t -\{eA- <LB)/t\ for (a) Ua/I = 1.0, and (b) (7aA = 5.0. 
Point of repulsion of flow lines on the Ub jt axis identifies a 
transition point between two types of phases. 



FIG. 5: RG flow diagram as projected on the UAjt — Uslt 
plane together with the phase boundaries (dashed lines) which 
seem to be the "lines of repulsion" of the flow lines. The "line 
of homogeneity" is shown by the dotted line. 



as well. Z/o found for two diflFerent types of sites need not 
be equal to each other, in general. 

To see the nature of the short range spin correlation 
we further compute the nearest neighbor (nn) spin corre- 
lation function l\S\zSiz)-, and the next nearest neighbor 
(nun) spin correlation function {SizS^z} in each type of 
cells, where Siz = rii^i — n^ i. Recursion relations are 
very much similar to that for Lq. 



V. RESULTS OF RENORMALIZATION GROUP 
CALCULATIONS 

We construct the phase diagram primarily by studying 
the RG flow pattern. We always start our iteration with 
= es = and UA/t,UB/t > 0. A transition point can 
be easily identified by looking into the RG flow diagram 
in the effective parameter space {Ub /t, \ {ea — eB)/t | } at a 
given value of UA/t. For each value of UA/t, there exists 
a point on the Us/t axis (with < Us/t < UA/t) which 
behaves like a "point of repulsion" between the flow lines 
(fig. 4). 

Starting from any point on its left the RG flow tends 
to go to the fixed point {0,oo}, indicating a charge or- 
dered phase, while any point on its right flows to {oo, 0} 
indicating a spin density wave. A plot of these transition 



points on the UA/t — UB/t plane shows the phase bound- 
aries, which can also be viewed as "lines of repulsion" 
of the flow diagram projected on the UA/t ~ UB/t plane 
(fig- 5). 

At this point it is interesting to find out the energy 
gap in the spectrum, A jig, from the present RG calcula- 
tion. For the present case we find that the Hamiltonian 
always flows under the RG iterations to flxed points cor- 
responding to the "atomic limit" {t — > 0) of (5) with 
parameters and C/j^°°^ = = U^°°\ where 

the superscript (oo) refers to the converged values of the 
corresponding parameters in (5). Therefore, the gap can 

be calculated from Arc = |(|e^°°^ - e^°°^| - C/(°°))|. We 
find that at the transition points the energy gap in the 
spectrum vanishes. 

In fig. 6 we plot the energy gap Arc as a function 
of Us/t for different values of UA/t. It is interesting to 
note that the zeroes of Ana occur at values of Ua and U b 
which are very close to the corresponding values obtained 
at zeroes of Ahfa- 

Therefore, the phase diagram obtained in the RG cal- 
culation agrees fairly well with that obtained in the HFA 
calculation (fig. 7). Departures are appreciable only at 
the strong coupling limit; it can be easily understood 
because in this limit the error due to the truncation of 
basis becomes most serious in the present RG scheme for 
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FIG. 6: Plot of the energy gap Arc, as found from the RG 
calculation, against Us/t- The dotted, solid, and the dashed 
lines correspond to Ua = 0.4, 1.2, and 2.5 respectively (scale 
of energy is chosen by setting t — 1.0). 

a model having two types of cells [l^ . The energy scales 
used for truncation of basis in the two types of cells now 
become appreciably different; the RG result, therefore, 
may not be highly reliable in this sector. However, the 
agreement of the RG and the HFA results from the weak- 
to the intermediate coupling regime really indicates that 
these results are very much reliable. 

In fig. 8 we plot the local moments at two types of 
sites as functions of C/s/t for different values of C/^/t. It 
turns out that Lq is always higher at sites with larger 
values of Coulomb correlation. This is in agreement with 
a previous observation in Q. Lq a.t A type site equals 
that at B type site only at the point of homogeneity i.e. 
at Ua — Ub- It is rather interesting to note that the 
values of Lq is slightly larger than | even in the param- 
eter space where, according to the RG flow pattern (and 
also from the HFA), the system develops a CDW instabil- 
ity. This shows that this correlation-driven CDW phase 
is dominated by short length-scale fluctuations that are 
suppressed in a CDW phase generated by a periodic mod- 
ulation in the site potentials alone p^ . 

We also calculate the nn and nnn spin-spin correlation 
functions with reference to ABA and BAB type cells. 
Plotted as functions of UB/t these (fig. 9) reveal some 
interesting points. The negative value of {SizS2z) in- 
dicates a nearest neighbor antiferromagnetic alignment, 
which is increasing from the CDW to the SDW regime, 
as one should expect on physical grounds. That this cor- 
relation is not zero even within the CDW phase is re- 
ally a reflection of the short-range fluctuations that we 
have just mentioned. However, the nn correlation has 
slightly different values in ABA and BAB cells. This dif- 
ference, which is much more pronounced at higher values 
of \Ua — Ub\, is possibly related to the finite size effect of 
the RG. This effect is suppressed at weak coupling and 
near the "homogeneous" point in the SDW region. On 
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FIG. 7: Phase diagram obtained from the RG (solid line) 
superimposed on the HFA phase diagram (dashed line) in the 
UA/t — Us/t plane. The "line of homogeneity" is shown by 
the dotted line. 



the other hand, the nnn correlation has opposite signs 
for the ABA and the BAB type cells in a region cor- 
responding to the charge density instability. This indi- 
cates a frustration that suppresses the large distance anti- 
ferromagnetic correlation. This frustration persists even 
within the SDW phase as it is evident from the widely 
different values of the nnn correlation in ABA and BAB 
type cells. Such an effect has already been anticipated 
in a previous work Q. It is to be noted here that this 
effect gradually reduces qsUb ^ Ua- In this respect the 
"SDW phase" of the half-filled AHM is behaving in a dif- 
ferent way than the SDW instability found in a half-filled 
Hubbard model (the "homogeneous limit" of the present 
model (1)). 



VI. STUDY OF MAGNETOCONDUCTANCE 

Now we study the nature of the ground state of the su- 
perlattice structure in presence of a magnetic field. The 
magnetic field penetrating the ring will interact with the 
moments of the electrons. There will be an additive Zee- 
man term of the form ii^H' in the Hamiltonian (1), where 
fia is the moment of the spin a and H' is the penetrating 
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FIG. 8; Plot of the local moments Lo as functions of Us/t 
for (a) Ua = 0.4, (b) Ua = 1.2, and (c) Ua = 2.5 (scale of 
energy is chosen by setting t = 1.0). The solid line shows the 
data for the ABA block, and the dashed line shows the data 
for the BAB block. The curves cross over each other at the 
"point of homogeneity" . 

magnetic field. The Hamiltonian now becomes 

z,cr z,cr 

+C/yi ^ni^tn^^x + ^«i,T"-i,i • (17) 

We follow the same HFA to decouple the Hamiltonian. 
Since we have found regions of the parameter space where 
the HFA results and the RG results match very well, we 
can rely on the HFA results in these regions. The Hamil- 
tonian corresponding to the up (down) spin electrons will 
generate the up (down) spin band. These two bands are 
not degenerate because the spin reversal symmetry is now 
broken. Therefore, at half-filling, the number of up and 
down spin electrons will be different and a net moment 
will be generated. 

For non-zero values of Ua{> 0) and Ub = 0, the system 
is charge ordered at half filling in absence of a magnetic 
field. This is an insulating phase with a gap between the 
upper and lower bands; at half filling the lower bands 
for both spin species are totally filled. When a small 
magnetic field is turned on, the up and the down spin 
bands shift in the opposite directions in energy scale. At 
a sufficiently large value of the magnetic field the up- 



FIG. 9: Plot of the spin-spin correlation functions as functions 
of UB/t {UA/t — 2.5) for (a) nn correlation, and (b) nun 
correlation. The solid line shows the data for the ABA block, 
and the dashed line shows the data for the BAB block. 



per band of the up spin electrons tend to get occupied 
at the cost of depopulating the lower band of the down 
spin electrons. Thus the ground state contains unequal 
number of up and down spins and the system, now hav- 
ing two partially filled bands, becomes conducting. If 
we go on increasing the magnetic field beyond a certain 
value, the said up spin band will be completely filled 
while the down spin band will be completely empty. The 
system is now spin polarized, leading again to an insulat- 
ing phase. This can be easily checked by calculating the 
Drude weight [T^ which measures the dc conductivity 
of the chain. Similar analysis can be made for Ub 
as well. Of course, in this case there arises a possibility 
of transition to a conducting phase from an SDW under 
the application of the field H'. We have calculated the 
Drude weight both for Ub — and Ub = 2. To calculate 
the Drude weight, a vanishingly small magnetic Hux $ 
(in units of flux quantum (f>o — hc/e) is introduced. The 
ring encloses this flux, but the flux does not penetrate 
the ring. Now the hopping term is modified by a phase 
factor. The Hamiltonian becomes 

1. a 

2, (7 
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FIG. 10: Plot of the Drude weight D vs. the field H' for (a) 
C/s = 0. and {h) Ub ~ 2.0. The dotted, dashed and the solid 
lines correspond to Ua = 1.0, 2.0, and 4.0 respectively in (a), 
and to Ua ~ 3.0, 2.0, and 0.9 respectively in (b) (scale of 
energy is chosen by setting t — 1.0). 



FIG. 11: Plot of the Drude weight D vs. system size A'' for 
(a) Ua = 1.0 and Ub = 0, and (b) Ua ^ Ub = 2.0 (scale of 
energy is chosen by setting t = 1.0). Cases with H' — (solid 
line) and H' > H'^ (dashed line) are shown where H'^ is the 
value of H' for which the system enters a conducting phase 
from an insulating one. 



+UA^n^,^ni^l + UB^ni^^n^^l (18) 

where (j) = ^/N, N being the number of sites in the ring. 
The Drude weight is calculated from the formula [13, US] 

where £'($) is the ground state energy of (18) calculated 
within the HFA. 

We have plotted the Drude weight D against the pene- 
trating magnetic field H' (fig. 10) with N = 300 for fixed 
set of values of C/4 and Ub- The curves clearly show the 
transition between conducting and insulating phases. For 
lower and higher values of H' the Drude weight is zero 
i.e. the system is insulating. But for moderate values of 
H' , (depending on the value of on-site Coulomb repulsion 
energy) the Drude weight is quite large, implying a con- 
ducting phase. For higher values of Ua the conducting 
region becomes narrower. For different combinations of 
Ua and Ub values, we have plotted (fig. 11) the Drude 
weight D against system size N. For the values of the 
penetrating magnetic field H' , for which the upper band 
of the up-spins starts to depopulate and the lower down- 
spin band gets partially populated, the system goes to 



a conducting phase. In fact, the Drude weight remains 
unaltered as we increase the system size. 

For H' — 0, the present system is found to be either in 
a CDW or in an SDW phase (fig. 1) which is insulating. 
Therefore, the conductivity of such a system of macro- 
scopic size must be vanishing. Thus in this region, the 
Drude weight falls sharply with increasing system size. 
Similar things happen for small values of H' for which 
the system is yet to pass on to the conducting phase. To 
identify the nature of the insulating phases directly, we 
study the spin and charge density parameters, s and c 
respectively, together with the magnetization 

m = -(ni3,T - + nA,^ ~ nA,i) . 

Fig. 12 shows the plots of c, s, and m as functions of 
H' for different values of Ua and Ub- For Ua >> Ub 
we find that the increase in H' drives the system from 
a charge ordered phase (large c) to a phase where c, s, 
and m are nearly comparable; this indicates the exis- 
tence of a conducting state due to the absence of any 
kind of long range ordering. We have already noted in 
fig. 10 that in this region D is non-zero which is consis- 
tent with the above observation. Further increase in H' 
would completely depopulate the down-spin bands and 
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FIG. 12: Plot of the order parameters c, s, and m against 
the field H' . The solid line shows magnetization, the dotted 
line shows spin order parameter, and the dashed line shows 
charge order parameter for (a) Ub ~ 0.6 and (b) Ub ~ 1.5 
while Ua ~ 2.5 (scale of energy is chosen by setting t = 1.0) 
and TV = 100. 



consequently the up-spin bands would be filled up. This 
induces a sharp rise in m (together with a vanishing D) 
indicating a transition into an insulating spin polarized 
phase. For Ub ^ Ua{^ t), however, the initial transition 
is from a spin ordered phase to a metallic one. 

Thus we can construct the phase diagram (fig. 13) 
of the AHM in presence of a Zeeman field which shows 
the possibility of having an enhanced magnetoconduc- 
tance for moderate values of the repulsive interactions. 
In case of very large Coulomb interaction, however, no 
conducting phase appears, and a direct transition from 
an insulating SDW phase to a spin polarized phase is 
observed. 



VII. CONCLUSION 

Summarizing, we have investigated the one dimen- 
sional half-filled alternating Hubbard super lattice struc- 
ture at zero temperature using the Hartree-Fock approxi- 
mation and a real space renormalization group technique. 
Although both the methods are approximate ones, they 
are complementary to one other. The agreement between 
the phase diagrams obtained from these two different 
methods shows that the results are rather reliable. We 



FIG. 13: The HFA phase diagram of the AHM in presence of 
a Zeeman field H' for iV = 100 and for (a) Ub /t = 0, and (b) 
Ub/I = 2 



obtained two types of phases, dominated by charge and 
spin density wave instabilities respectively, depending on 
the Ua and Ub values. The Ua — Ub phase diagram 
shows two transition lines, indicating a centrally located 
antiferromagnetic region and two charge density wave re- 
gions near the axes. 

The system is insulating. It may be noted at this point 
that the possibility of such transitions was not explored 
in the previous works on this model; most of these works 
concentrated only on a limited region of the parameter 
space, e.g. Ua > and Ub — 0, and also there were 
some severe restrictions on the system size. Apart from 
the phase transitions, it is interesting to note that there 
appears a spatial modulation of the local moments dic- 
tated by the inhomogeneity of the correlation parame- 
ters. The RG results also indicate that the ordering is 
not full grown as compared to the homogeneous limit be- 
cause of underlying frustration of spin-spin correlation. 
These results are in agreement with a previous finding 
4]. Also we have studied the ground state properties, 
including the Drude weight, of the superlatticc in pres- 
ence of a magnetic field {H'). The Ua — H' (for fixed 
Ub) phase diagram shows that the system becomes con- 
ducting for an intermediate range of values of H' . In the 
phase diagram, we obtained insulating phases for lower 
and higher values of H' . The width of the conducting re- 
gion depends on the values of on-site Coulomb repulsion 
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energies of the super lattice structure. It seems inter- Acknowledgments 
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